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Abstract. In this paper HIPPARCOS astronietric and 
kinematic data are used to calibrate both infrared lumi- 
nosities and kinematical parameters of Long Period Vari- 
able stars (LPVs). Individual absolute K and IRAS 12 
and 25 luminosities of 800 LPVs are determined and made 
available in electronic form. 

The estimated mean kinematics is analyzed in terms of 
galactic populations. LPVs are found to belong to galac- 
tic populations ranging from the thin disk to the extended 
disk. An age range and a lower limit of the initial mass is 
given for stars of each population. A difference of l.Smag 
in K for the upper limit of the Asymptotic Giant Branch 
is found between the disk and old disk galactic popula- 
tions, confirming its dependence on the mass in the main 
sequence. 

LPVs with a thin envelope are distinguished using the 
estimated mean IRAS luminosities. The level of attraction 
(in the classification sense) of each group for the usual 
classifying parameters of LPVs (variability and spectral 
types) is examined. 

Key words: Stars: absolute magnitude ~ variable stars - 
kinematics - galactic populations - AGB 



1. Introduction 

Long period variables (LPV) form an important class of 
red giant stars. They show more or less regular photomet- 
ric variability with amplitudes reaching up to 8 magni- 
tudes and periods up to 600 days. They traditionally com- 
prise Miras, semi-regular (SR) and irregular (L) variables 
according to the amplitude and the regularity of their vi- 
sual light curves. They are known to be either 0-rich or 
C-rich, and comprise thus M, S and C stars. More recently 
OH-IR sources have been found from infrared and radio 
observations showing that they belong to the LPV popu- 
lation with periods up to 2000 days. Those sources emit in 
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the infrared and radio wavelengths, and are not associated 
with any detectable counterpart in optical wavelengths. 

The brightest LPVs are luminous enough to be ob- 
served at long distances, providing information on the host 
galaxy, like the Magellanic Clouds (see Van Loon et al., 
1999 as an example). While the ranges of masses and ages 
of LPVs are still the subject of discussion, it is generally 
accepted that they are large, and are therefore considered 
as very good tracers of galactic history. 

The determination of the characteristics of individual 
LPVs is usually a delicate task due to the complexity of 
the dynamic and chemical phenomena to be considered. 
A statistical study using all available data of a large 
sample of LPVs is often needed. A rough example of 
such an approach could be the relation between the mean 
visual light curves and the infrared colors of C and 0-rich 
LPVs already presented in Mennessier et al. (1997a). 
In this paper, HIPPARCOS astrometric data and the 
available multi-wavelength (K, IRAS 12 and 25) infrared 
photometric measurements allow us to calibrate multi- 
wavelength luminosities and to discriminate between 
different galactic populations - and thus different ranges 
of initial masses {Mms) - among the LPVs according to 
their kinematical properties. In a second step, individual 
K and IRAS absolute magnitudes are estimated for all 
the 800 considered LPVs using a powerful statistical 
estimator. 

Our sample of LPV stars and the data used are de- 
scribed in Sect. ^. The statistical method specifically de- 
veloped for the study of HIPPARCOS samples is summa- 
rized in Sect. ||. Sect. ^ presents the discriminated groups 
of LPVs resulting from our statistical analysis, while Sect. 
H analyzes the results derived for individual stars. Finally, 
Sect. I reviews the crossed properties derived from the 
analysis at different wavelengths. 



2. The sample of LPV stars 
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2.1. Data 



In order to benefit from the accurate astrometric data 
made available by the HIPPARCOS satellite, we use in 
our study the sample of all LPV stars observed on this 
mission, i.e. the LPVs brighter than 12.5 mag in V during 
more than 80% of their variability cycle. The sample is 
composed of about 900 stars which are either of type M 
(0-rich), C (C-rich) or S (0/C <, 1). They include Mira, 
SR (of both type a and b) and L variables. 

Astrometric data is taken exclusively from the HIP- 
PARCOS Catalogue (Perryman et al. 1997) to provide a ; 
homogeneous data set. Radial velocities are taken from the 
HIPPARCOS Input Catalogue (HIC; Turon et al. 1992). 

Photometric data are gathered from various sources. 
V magnitudes (mv) are taken from the HIC. They corre- 
spond to the magnitudes given in the General Catalogue of 
Variable Stars (GCVS; Kholopov et al. 1985), corrected as 
described in the HIC volumes to obtain mean magnitudes 
at the maxima of light. K magnitudes {mx) are taken 
from the Catalogue of Infrared Observations (Gezari et al. 
1996), and include the large set of JHKL measurements of 
LPVs by Catchpole et al. (1979) and the measurements by 
Fouque et al. (1992), Guglielmo et al. (1993), Groenewe- 
gen et al. (1993), Whitelock et al. (1994), Fluks et al. 
(1994), Kerschbaum & Hron (1994), Kerschbaum (1995) 
and Kerschbaum et al. (1996). Infrared magnitudes are de- 
rived from the F12 and i^25 fluxes measured at 12 and 25 
micrometers respectively by the infrared astronomy satel- 
lite (IRAS). We use 
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mi2 = 3.63 - 2.5 x logFi2 
and 

77125 = 2.07-2.5 X logF25, 



(1) 



(2) 



as given in the IRAS-PSC catalog (vol.1. Sect. VI.C.2). 
One should note that not all authors use this definition for 
the infrared magnitudes 77712 and 77125. The color index 
W25 — mi2 used by Van der Veen and Habing (1988), for 
instance, is higher than the one deduced from Eqs. |^ and 
^ by 1.56mag. 

Among the ~900 stars of our sample, the number of 
stars for which V, K and IRAS infrared magnitudes are 
available amounts to 882, 652 and 793, respectively, with 
608 stars having both K and IRAS magnitudes. 

Finally, variability and spectral types are taken from 
the GCVS. 

2.2. Selection effects 

The main selection bias in our sample comes from the HIP- 



PARCOS magnitude limit V < 12.5 mag (see Sect. 2.1) 



This selection is well determined and thus easy to take 
into account in the statistical analysis. 

The characteristics of LPVs cause another bias re- 
lated to the magnitude limit of the sample. LPV stars are 
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Fig. 1. Distribution of the GCVS LPVs according to 
apparent visual magnitude at maximum luminosity and 
IRAS color index. Stars observed by HIPPARCOS are in- 
dicated 



evolved red giants, often characterized by the formation 
of dust around them. The presence of a dusty circum- 
stellar envelope affects the stellar spectrum by reducing 
their visible light. As a result, obscured LPVs are under- 
represented in our sample, because the HIPPARCOS se- 
lection was done on the basis of the visible magnitude. 
The importance of this bias can be estimated by com- 
paring the number of stars included in the HIC with the 
number recorded in the GCVS. This comparison is shown 
in Fig. 1^ as a function of the V magnitude 77iy and the 
color index 77125 ~ ™i2j where we consider all LPVs of 
the GCVS for which either the visual (777y) or the pho- 
tographic (77ip) magnitude at maximum is given, and as- 
suming mp — 77iy = 1.8 as the mean value for LPVs. All 
stars from the HIC, represented by filled circles in Fig. ^ 
are found to have 777y > 12.5, as expected (the very few 
exceptions being most probably due to the fact that the 
assumed mp — mv = 1.8 relation does not apply to them). 
Fig. also shows that the number of stars included in 
the HIC (relative to the number of stars recorded in the 
GCVS, represented by filled and open circles in Fig. |l]) 
decreases with increasing circumstellar envelope thickness 
(i.e. decreasing 7772^— 7(1,12 index). This bias is further dis- 
cussed in Sect. 5.2.3 . 

It must be noted that the GCVS itself is, of course. 



not exhaustive, and is certainly biased at the expense of 
the reddest stars. OH-IR stars, for example, are not well 
represented in the GCVS sample. For these reasons, a sta- 
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tistical method which can take into account aU these biases 
is necessary for our analysis. This method is described in 
the next section. 

3. The statistical method 

In this paper the Luri-Mennessier (LM) statistical 
method, described in Luri et al. (1996), is used to analyze 
the sample of LPV stars. The method has been specifi- 
cally designed to exploit the HIPPARCOS data and thus 
is suitable for our purposes. This method has already given 
fruitful results, in particular for Barium stars (Mennessier 
et al., 1997c), Ap-Bp stars (Gomez et al., 1998) and the 
LMC distance modulus (Luri et al., 1998). 

The use of appropriate statistical methods for the ex- 
ploitation of the HIPPARCOS astrometric data is crucial 
in order to obtain correct results. Otherwise the values ob- 
tained can be affected by strong biases and the precision 
of the data will not be fully used. A discussion on the cor- 
rect use of HIPPARCOS data and recommendations on 
analysis techniques can be found in Brown et al. (1997). 

The LM method used in this paper is especially ap- 
propriate for use with the HIPPARCOS data. We refer 
to Luri et al. (1996) for a detailed description, and only 
briefly summarize here some of its main features. 

First of all, the stellar population from which the sam- 
ple is extracted is assumed to be composed of several dis- 
tinct groups. These groups can differ in kinematics, lumi- 
nosity or spatial distribution and its number is a priori not 
known. Therefore, using a sample extracted from this base 
population and taking into account the selection criteria 
used to create it, the LM method: 

— determines the number of significant discriminating 
groups and the percentage of each of them in the pop- 
ulation; 

— produces, for each group, unbiased estimates of the 
parameters of the model used to describe it, i.e.: 

— a Schwarzschild ellipsoid velocity distribution char- 
acterized by {UQ,Vo,Wo,au,crv,crw)- Uo,Vq and 
Wo point towards the galactic center, the galactic 
rotation, and the north galactic pole, respectively; 

— an exponential distribution of the number of stars 
in the direction perpendicular to the galactic plane, 
with a scale height Zq; 

— a normal distribution of the absolute magnitude 
in the bandpass of the used observed magnitudes 
characterized by the parameters Mq and aM', 

The results of this first step for our LPVs sample are 
given in Sect. 

The minimum input data needed by the LM method 
are the measured positions, proper motions and apparent 
magnitudes of the stars, but it can also use the paral- 
laxes and radial velocities if available. The method takes 



into account the selection effects of the sample, the obser- 
vational errors, the galactic rotation and the interstellar 
absorption. 

In a second step, once the groups are identified and 
parametrized, the method: 

— computes for each star of the sample, the a posteriori 
probability that the star belongs to a given group; 

— uses a Bayesian rule to assign each star to a group; 

— uses a statistical estimator to obtain estimations of 
individual distances and luminosities for each star. 

This second step is presented in Sect. || where the rep- 
resentativity of our sample with respect to the kinematic 
and photometric properties of the population is also dis- 
cussed. 

4. CaUbrations of the LPVs population 

The LM method was applied four times, as described in 
Sect. H, once for each photometric bandpass: V - results 
already presented in Mennessier et al. (1997b) -, K, 12 and 
25 - in the present paper. In principle, one could assign 
a joint luminosity distribution to two or more bandpass 
magnitudes simultaneously, and the LM method would 
separate the sample into stellar groups consistent with 
all the photometric measurements together. This option, 
however, requires a perfectly well known relationship be- 
tween the different magnitudes in order to define a joint 
distribution function as realistic as possible for all the 
band passes. The correlation between the near-infrared 
(K) and IRAS infrared properties presently cannot be well 
modeled and very likely has a non-unique form depend- 
ing on the stellar and circumstellar evolutive stage along 
the AGB. Thus, we decided not to couple the photometric 
band passes and to calibrate each luminosity separately. 
Furthermore, bandpasses are related to different physical 
processes and can provide separate interesting informa- 
tion: V is greatly affected by absorption molecular lines, K 
refiects the stellar emission, and IRAS bandpasses depend 
on the nature and density of grains in the circumstellar 
envelope. 

The LM method is simultaneously sensitive to kine- 
matics and luminosity and thus the number of significant 
discriminating groups depends on both these characteris- 
tics and is not necessarily the same for the different band- 
pass analyses. Furthermore, the samples used are not the 
same, and this can also affect the number of discriminated 
groups. 

Six distinct groups are identified in the V magnitudes, 
three in K and four in each of the two IRAS magnitudes. 
Those are successively analyzed in terms of the classical 
galactic populations. Although the number of groups is 
found to be different for each analysis, the groups present 
similarities in their kinematical composition and with re- 
spect to the galactic populations (see sect. ||). 
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Table 1. V calibration: estimated parameters of the dif- 
ferent groups and percentage of the sample into each of 
them. 



Group 


BD 


D 


ODl 


0D2 


TD 


ED 


Mv 


-3.6 


-1.0 


-1.2 


-0.2 


-1.2 


-2.8 


a My 


1.4 


0.8 


0.2 


1.0 


0.5 


1.2 


Uo 


-10 


-6 


-44 


-1 


-34 


-61 


Vo 


-11 


-6 


-35 


-21 


-84 


-235 


Wo 


-13 


-6 


-6 


-10 


-19 


-20 


cruo 


13 


24 


28 


37 


77 


188 




14 


14 


25 


23 


29 


126 


CWo 


9 


9 


22 


23 


65 


72 


Zo 


104 


126 


217 


249 


409 


1227 


% 


8 


25 


13 


44 


8 


2 



Table 2. K calibration: estimated parameters of the dif- 
ferent groups and percentage of the LPVs population in 
each of them. 





Group D 


Group OD 


Group ED 




est. 


a 


est. 


a 


est. 


a 




-6.1 


0.4 


-6.0 


0.7 


-5.3 


0.8 


(TK 


1.1 


0.3 


0.7 


0.4 


1.4 


0.5 


Uo 


-7 


10 


-17 


27 


-21 


14 


(^u 


29 


9 


45 


16 


111 


11 


Vo 


-12 


8 


-36 


25 


-123 


12 


av 


16 


5 


27 


11 


69 


18 


Wo 


-9 


6 


-6 


9 


-20 


19 


aw 


12 


3 


26 


13 


90 


25 


Zo 


184 


44 


268 


85 


782 


313 


% 


60 




35 




5 





4.1. The V band 

An analysis of the six groups identified in the V band 
has been presented in Mennessier et al. (1997b). In order 
to compare these results with the ones obtained for in- 
frared calibrations, the main results are summarized here. 
Table reviews the estimated mean parameters for the 
analysis corresponding to the V luminosity at the phase 
of maximum light. The LPVs are found to belong to all 
galactic populations from disk to very extended disk. We 
wish to emphasize three points: 

— LPVs belonging to the bright disk population (BD) 
have a mean luminosity of My=-3.6 and Zo=104pc, 
while the corresponding values for the disk population 
(D) are My=-1.0 and Zo=126pc. Among LPVs be- 
longing to BD population there are probably stars at 
the upper limit of the AGB or even at the first step of 
the post-AGB state. 

— the main group (44%) has an estimated scale height 
of Zo=249pc. Its kinematics is typical of the old disk 
population. 

— a few LPVs (less than 2%) belong to the extreme ex- 
tended disk (ED), with Zo greater than 1200pc and a 
very large Schwarzschild ellipsoid velocity distribution. 
This most probably indicates a birth of those stars 
early in the evolution of our Galaxy. Those stars should 
thus be a metal-deficient disk population (see |6.2|) . 
This is consistent with the bright luminosity found 
(My =-2. 8). 



4.2. The K hand 

Only three groups are identified in the K band. From their 
kinematics and spatial distribution, given in Table ^, they 
can be interpreted as the galactic disk (D), old disk (OD) 
and extended disk (ED) populations. They are simil ar t o 
the four main groups identified in the V band (Sect. |4.l|) , 



Table 3. 12 calibration: estimated parameters of the dif- 
ferent groups and percentage of the LPVs population in 
each of them 





Group D 


Group ODb 


Group ODf 


Group ED 




est. 


a 


est. 


a 


est. 


a 


est. 


a 


12o 


-6.4 


0.3 


-8.0 


0.4 


-6.4 


0.5 


-6.2 


1.0 




1.7 


0.1 


1.2 


0.2 


0.6 


0.4 


1.6 


0.2 


Uo 


-6 


9 


-12 


5 


-10 


9 


-30 


37 


(TU 


22 


6 


35 


8 


39 


6 


106 


45 


Vo 


-7 


8 


-26 


7 


-24 


8 


-97 


54 


av 


12 


5 


26 


11 


22 


6 


65 


64 


Wo 


-9 


4 


-9 


5 


-8 


5 


-2 


44 


aw 


9 


8 


21 


8 


21 


7 


75 


29 


Zo 


161 


55 


258 


56 


256 


79 


1065 


724 


% 


29 




32 




29 




10 





except that the disk and a part of the old disk population 
seem to be mixed. 

In a previous analysis of a sample restricted to O- 
rich Miras (Alvarez et al.,1997), only two groups were 
found. One corresponded to the extended disk popula- 
tion, with a percentage of 17, in agreement with our re- 
sult i.e. 18/142=13% of 0-rich Miras belonging to the ED 
group (see table 0) . The other group mixed disk and old 
disk populations. In the present paper, a more numerous 
sample allows a more refined separation of the kinematic 
populations. 

4.3. The IRAS bands 

The four LPV groups identified in the IRAS 12 and 25 
bands are given in Tables ^ and ^, respectively. They are 
similar to those identified in the K band (Table ||), except 
that the old disk group is further divided into "bright" 
(ODb) and "faint" (ODf) subgroups. Let us remember 
here that the method allows us to distinguish groups with 
similar mean kinematics but different luminosities (ODb 



M.O. Mennessier et al.: Photometric and kinematic properties of LPVs 



5 



Table 4. 25 calibration: estimated parameters of the dif- 
ferent groups and percentage of the LPVs population in 
each of them 





Group D 


Group ODb 


Group ODf 


Group ED 




est. 


a 


est. 


a 


est. 


a 


est. 


a 


25o 


-7.1 


0.5 


-8.6 


0.4 


-6.5 


0.3 


-6.8 


0.8 


0"25 


1.7 


0.1 


1.2 


0.2 


0.6 


0.4 


1.6 


0.5 


Uo 


-6 


6 


-10 


4 


-10 


7 


-39 


48 


cru 


21 


8 


36 


10 


38 


6 


111 


33 


Vo 


-6 


4 


-26 


7 


-22 


6 


-99 


63 


cry 


13 


4 


27 


9 


22 


5 


69 


23 


Wo 


-10 


4 


-9 


5 


-8 


4 


1 


42 


aw 


11 


7 


21 


8 


20 


4 


75 


37 


Zo 


158 


43 


277 


34 


270 


107 


1610 


1180 


% 


28 




32 




30 




10 





nitudes. We can deduce that the accuracies of our esti- 
mated individual luminosities are distributed according 
to a gaussian rule of standard error 0.3 mag. and 0.1 mag. 
respectively in K and IRAS bands. The lower accuracy 
in the IRAS bands is consistent with the fact that IRAS 
photometry is more homogeneous than the K photometry, 
and that the variability amplitude of LPVs is smaller in 
the IRAS bands than in K. As previously stated, the LM 
method has allowed us to take advantage of all the avail- 
able information, leading to better estimations of the indi- 
vidual absolute magnitudes. Furthermore, the LM method 
has provided at the same time the statistical distribution 
of these individual magnitudes (see Sect. ^ and The in- 
dividual estimates of K, 12 and 25 luminosities arc given 
in a table available in electronic form at the CDS f| and 
are included in the specialized ASTRID database 0. 



and ODf) or groups with a similar luminosity distribution 
but different kinematics (D and ED for instance). More- 
over, it is important to remark that D, ODb and ED have, 
on average, a similar color index 25 — 12 = —0.6 mag corre- 
sponding to a thick circumstellar envelope, while ODf has 
a mean index of 0.1 mag that suggests that the majority 
of the stars in this last group have thin envelopes. 

Let us finally point out that the kinematic parameters 
{Uo,Vo,Wo) associated with each of the four groups are 
very similar for both the 12 and 25 calibrations. 



5.2. Comparison of sample/population 

The LM method gives unbiased calibrations for the base 
population. It also gives individual kinematic and photo- 
metric estimates for each star of the sam ple. The distribu- 
tion of these individual estimates (Sect. pT| ) is, of course, 
biased by the sample selection criteria, contrary to the 
group characteristics derived in Sect. ^. A comparison of 
the statistical properties of the sample with the calibrated 
parameters for the population allows us to check the rep- 
resentativity or the bias of the sample with respect to the 
population. 



5. Individual estimates and properties of the 
sample 

5.1. Individual estimates 

Once the parameter estimation and group discrimination 
is completed, each star in our initial sample is a poste- 
riori attributed to one of the LPV groups identified in 
each bandpass, following the method described in Sect. ||. 
This allows us to estimate the most probable individual 
distance and absolute magnitude in each band according 
to the observed astrometric, kinematic and photometric 
data and attributed group. 

Due to the probabilistic nature of the Bayesian proce- 
dure, some misclassification is unavoidable. To check and 
improve individual star assignations in each wavelength, 
we compare the calculated color indices cal = — 
(obtained from the estimated individual absolute mag- 
nitudes deduced by the Bayesian assignations in the Ai 
and/or A2 wavelengths) with the observed color indices 
(obs = — mx^). 5% of the stars are re-assigned to 
groups reducing the differences cal — obs for their indices 
25-12 and/or K-12. 

Figure ^ shows the histograms of difference of cal — obs 
for the indices 25-12 and K-12 of all the stars in the sam- 
ple. These distributions are related to the errors of the in- 
dividual estimated luminosities and of the observed mag- 



5.2.1. Kinematic representativity 

Let us analyze the representativity of our sample stars 
with respect to the kinematics. The observed proper 
motions and radial velocities, together with the estimated 
individual distances, allow us to compute the three 
velocity components {U, V, W) and the distance Z above 
the galactic plane for each star in our sample. The 
mean kinematical properties of our sample derived from 
these individual estimates are shown in Table ^ for each 
group of the K and 12 bands. They are very similar 
to the parameters describing the groups (Tables || and 
^), showing that our sample is very representative of 
the LPVs population as far as the kinematics is concerned. 

This conclusion was expected since there is a priori no 
selection affecting (directly or indirectly) the kinematics 
of our sample and thus no bias is introduced in the kine- 
matics of the stars. We can also note that the proportion 
of the different groups in the sample is close to that found 
in the population. 



^ via anonymous ftp to cdsarc.u-strasbR.fr (130.79. 128.5) or 
via http://cdsweb.u-strasbs.fr/csi-bin/qcat? J/ A+ A/ 



tp, 

via pittp : / / astrid.graal. univ-montp2 . fr 
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Fig. 2. Histograms of the distributions of the differences of the observed color indices (obs) and the calculated (cal) 
arbitrary scale - from estimated intrinsic luminosities for initial and final discriminations of stars into the groups 



LPV's - LM calibrations 



LPV's - LM calibrations 



LPV's - LM calibrations 




n m n , ^^.^ 




LPV's - LM calibrations 



LPV's - LM calibrations 



LPV's - LM calibrations 



LPV's - LM calibrations 



D 


SO.0 


JDh 


SO.0 




.TDf 

L 




.ED 


















LPV s - LM calibrations 


LPV's - LM calibrations 


LPVs - LM calibrations 


LPV's - LM calibrations 


,-D 


SO.0 


.TDb 


SO.0 




.TDl 




.ED 


i.... 



















Fig. 3. Distributions of individual luminosities in K, 12 and 25 from top to bottom, of each group (D,OD ED, or 
D,ODb,ODf,ED from left to right) of the sample compared to the distribution of the calibrated luminosity (in units 
normalized to the surface of each histogram) for the same group 



5.2.2. Luminosity representativity 

Although no kinematical bias is introduced when se- 
lecting a sample with a cut in apparent magnitude, it 
is well known that a bias in luminosity is introduced. 
Figure ^ shows the histograms of the individual absolute 
magnitudes of the stars in our sample in each group of 



the K and IRAS bandpasses, together with the normal 
unbiased distributions estimated by the LM method 
for the population. These distributions are in units 
normalized to the surface of each histogram -and not to 
the number of population stars in each sample -, thus 
only the relative shapes and the magnitude shifts of both 
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histogram and unbiased distribution are relevant. The 
bias of our sample towards higher luminosities is very 
clear both in K and IRAS bands. In the IRAS bands, 
the under-representativity of faint stars in our sample is 
more pronounced for LPV stars in the disk group than in 
the other IRAS groups. This corresponds to the classical 
Malmquist bias (1936), increasing with increased ctm 
value. In short, the under-representation of faint stars in 
our sample is important for the K or IRAS faint stars 
and even more for the disk population, specially in the 
case of IRAS bandpasses. 

However, let us remark that the brightest stars in every 
group of the sample coincide with the brightest luminosity 
of the group base population. 

5.2.3. Envelope effects and representativity 

The luminosity sampling bias is not independent of the ex- 
istence, thickness and composition of a circumstellar enve- 
lope around LPVs. Figure ^, which shows the percentage 
of known LPVs measured by HIPPARCOS (LPVs:%HIP) 
as a function of the IRAS (25-12) color index, shows that 
the incompleteness depends on the IRAS color. This is not 
surprising because the thicker the envelope, the fainter the 
star in the visual wavelengths. 

This is confirmed if instead of using the known LPVs 
we use the IRAS sources with a (25-12) color index 
compatible with the LPVs values of this index. In doing 
so, we include stars in the LPV region not necessarily 



Fig. 4. Percentages of LPVs observed by HIPPARCOS 

(full circles) as a function of the 25-12 IRAS color index. 
They are compared to the percentages of stars observed 
by HIPPARCOS (+) and of known LPVs (*) among the 
sample of IRAS sources selected as probable LPVs from 
their IRAS color indices. 

classified as variables (IRAS sel.:%LPVs). The bias of 
the HIPPARCOS sample is more strongly dependent 
on the envelope thickness if we do the comparison with 
these selected IRAS sources. Thus the percentage of stars 
observed by HIPPARCOS (IRAS sel.:%HIP) strongly and 
abruptly increases up to 80 % for 25-12 decreasing to zero. 

Finally, figure ^ shows how much the sample of carbon- 
rich stars measured by HIPPARCOS (C stars:%HIP) does 
not represents either the percentage of the C-rich stars 
among the known LPVs (LPVs:%C stars) or the percent- 
age of stars known as LPVs measured by HIPPARCOS 
(LPVs:%HIP). Thus one should be careful about making 
any interpretation from the percentages of C-rich stars, as 
we will see in sect. |6.4| . 

6. Properties of crossed groups 

6.1. Crossing K and IRAS 

Due to the difficulties coming from the large uncertainties 
on TOy: large amplitudes and variations from one cycle 
to another, the V results will not be further used in our 
analysis and we will concentrate on the K and IRAS 
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Table | gives the number of stars in our sample which 
are assigned to every crossed group G(G') i.e. to group 
G in K and G' in IRAS. In Sect.|5.2.l| our LPVs sample 



100.0 



■ LPV's:°/, 
^ LPV's:°/. 
3Cstars;%HIF 



60.0 




Fig. 5. Percentages of C stars observed by HIPPARCOS 

(empty circles) and among known LPVs (empty squares) 
as a function of the 25-12 IRAS color index compared with 
the percentages of LPVs observed by HIPPARCOS (fuU 
circles). 



results. 

The luminosity estimations in the K and IRAS bands 
complement each other in the sense that, in general, K 
fluxes characterize stellar properties while IRAS fluxes 
provide information on the circumstellar envelope. Thus, 
the most physically interesting results are obviously ob- 
tained by simultaneously considering K and IRAS lumi- 
nosities. We have already seen that it is very difficult to 
calibrate these luminosities at the same time due to the 
incomplete knowledge and non-uniqueness of the relation 
between the different magnitudes (Sect. ^). Another pos- 
sibility is to do a crossing of the groups from both the 
K and IRAS calibrations i.e. to examine the properties of 
the stars belonging to the same group in K and IRAS. 

The first remarkable result concerns the number of 
crossed groups: only 7 are not empty while 12 could, a 
priori, be expected. Interestingly, there is no mixture of 
the extended disk group (in either wavelength) with any 
other group, except two stars (0-rich SRa: RW Eri and 
0-rich Mira: SV And) , which is compatible with the statis- 
tical classification errors. This is a nice confirmation of the 
power of the LM method to extract consistently distinct 
groups in biased samples of a given stellar population. 



is shown to be representative of the LPVs population 
as far as the kinematics is concerned. Thus the mean 
kinematics of the stars belonging to a crossing group 
K(IRAS) may be considered as representative of the mean 
kinematics of the LPVs population belonging to this group. 

Obviously such a co nsider ation does not apply to 
the luminosities (see Sect. 5. 2. 2 ). The assigned groups are 
given in annex A (electronic table). 

6.2. Ages and initial masses 

Table || gives the values of the axes of the velocity 
ellipsoids and the scale height of each of the 7 crossed K 
and IRAS groups. Given that our sample is representative 
of the population in terms of kinematics, as already seen 
in Sect. 5.2.1, we can use the kinematical values of table 



^ as representative in terms of galactic populations. 

The relation between the mean kinematics of a galactic 
population and its age allows us to estimate the range of 
ages of the groups. Furthermore, classical statistical stud- 
ies of stars known to belong to different galactic popu- 
lations and of different metallicity abundances allow us 
to add an estimate of the range of metallicity. By com- 
paring the values in table ^ with the results on kinematics 
and metallicity of the galactic populations by Mihalas and 
Binney (1981) and by Stromgren (1987), we can deduce: 

— D(D) and OD(D) populations are 1-4 10^ yr old with a 
solar metallicity. Part of the stars in D(D) are classified 
as belonging to the bright disk population (BD) by 
the V analysis, being younger and with a small over- 
abundance. 

— The age of D(ODf) and D(ODb) populations can be 
estimated to be in the range 4-8 10^ yr, with a metal- 
licity from the solar one to [Fe/H] around -0.4 i.e. Z 
between 0.006 and 0.02. 

— OD(ODf) and OD(ODb) populations are composed of 
stars older than 8 10^ yr, up to more than 10^° yr, 
with [Fe/H] from the solar value to -0.7 i.e. Z between 
0.004 and 0.02 

— Stars classified as belonging to ED are probably very 
old and deficient with [Fe/H] between -0.7 and -1.5 i.e. 
Z of the order of (0.001,0.004). 

Moreover, we can estimate initial masses from evolu- 
tionary tracks. From Binney and Merrfield (1998) we can 
estimate a lower limit of the initial mass of a given age 
star that has reached the AGB. Thus, values of 2, 1.4, 
1.15, 1 M.Q can be deduced as lower limits of TWms of the 
stars of solar metallicity of respectively 1, 4, 8, 12 10^ yr. 
This agrees with the results on the ages at the top of the 
early- AGB (Charbonnel et al.,1996). 
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Table 6. Mean kinematical parameters for the crossing K(IRAS) groups computed from individual velocities and 
positions. They can be considered as representative for the LPVs population (see Sect. 6.1 





Diskl 


Disk2 


Old Disk 


Ext. Disk 




D(D) OD(D) 


D(ODf) D(ODb) 


OD(ODf) OD(ODb) 


ED(ED) 


nb of stars 


141 21 


103 90 


81 113 


36 


Vo 


-6 -7 


-18 -19 


-34 -32 


-123 


cru 


23 28 


35 36 


42 40 


114 


(TV 


13 23 


20 18 


24 26 


63 


aw 


11 11 


34 16 


21 29 


80 


Zo 


166 191 


208 231 


160 310 


620 


age range 


1-4 10" yr 


4-8 10" yr 


8-gtlO 10" yr 




lower mass limit 


2-1.4 Mq 


1.4-1.15 A^Q 


1.15-ltl Mq 





All these results are in the same ranges as the ones 
given by Jura and Kleinmann (1992), but our classification 
is more refined because it is not based on the spectral 
types and periods which are now known to not really be 
discriminative parameters for LPVs. 

In the rest of this paper the stars belonging to 
D(D) and OD(D), D(ODf) and D(ODb), OD(ODf) and 
OD(ODb), ED(ED) will be called diskl, disk2, old disk 
and extended disk LPVs respectively (see tables ^ and ^ . 

6.3. Evolutionary tracks 

Figure ^ shows the K magnitude as a function of the V- 
K color index for each of the diskl, disk2, old disk and 
extended disk groups. For comparison, evolutionary AGB 
model predictions are also shown for three different masses 
(1.5, 2.5 and 4Mq) at solar metallicity, and at three dif- 
ferent metallicities (Z=0.004, 0.008 and 0.02) for 2.5 Mq 
stars. These models have been computed at Geneva, and 
are described in Mowlavi (1999) and Mowlavi & Meynet 
(2000). The conversion between model variables (effective 
temperature Teff and luminosity L) to observable quan- 
tities {V — K and Mk) was done by using the transfor- 
mations given by Ridgway et al. (1980). 

Several uncertainties affect both model predictions and 
the color transformation relations for AGB stars (which 
are characterized by peculiar chemical compositions as a 
result of dredge-up episodes). They also affect the deter- 



mination of the stellar V magnitudes as noted in Sect. 4.1 



Thus, the comparison between the evolutionary tracks and 
the distributions of our sample stars in each group shown 
in Fig. ^ can only provide qualitative results. 

Evolutionary tracks show that, at a given metallicity, 
stellar luminosities increase with initial stellar mass (at a 
given V-K). We thus conclude, at least qualitatively and 
due to the solar or slightly deficient abundance of the ma- 
jority of HIPPARCOS stars, that our sample stars have 



lower mass limits which respectively decrease as we con- 
sider diskl, disk2 and old disk groups. This is in agreement 
with the conclusions drawn in Sect. |6.2| . Stars of the ex- 
tended disk group, on the other hand, are compatible with 
lower metallicities, given their higher K luminosities. 

We note that the evolutionary tracks cannot, even if 
freed from any uncertainty, attribute a single (M,Z) set of 
parameters to a star because of the degeneracy of those 
two parameters. A higher luminosity in K at a given V- 
K could be either attributed to a higher initial mass or 
lower metallicity. Kinematics can help in distinguishing 
such ambiguous cases. 

Finally, we can comment on the smaller V-K values for 
the bright disk LPVs. This confirms the strong circumstel- 
lar absorption in V for these massive stars. 

6.4- Variability and spectral types 

The composition of each crossed group K(IRAS) with re- 
spect to usual classifications of LPVs (variability and spec- 
tral types) is given by the contingency table of both assig- 
nations (Table 0) . The associated attraction-repulsion in- 
dices (Tenenhaus, 1994) - ratio of the observed frequency 
to the theoretical frequency in the case of independence 
of both modalities - are more significant in characterizing 
the correspondence analysis of types and of groups. These 
indices are given in Table attract or repel each other if 
the attraction-repulsion index is larger or smaller than 1 
respectively. 

^From Table ^ we can deduce that: 

— The groups corresponding to initially less massive stars 
are essentially attractive for 0-rich LPVs and the ones 
corresponding to initially massive stars are essentially 
attractive for G-rich LPVs. 

— The strong attraction of C-rich stars by groups cor- 
responding to more massive initial stars is clear. It 
agrees with the results already given by V calibration 
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Fig. 6. Theoretical evolutionary AGB tracks for stars of 1.5, 2.5 and 4 Mq with solar metallicity and for deficient 
(Z=0.008 and Z=0.004) stars of 2.5 compared to distribution of the individual estimated K luminosities as a 

function of the color index V-K according to the assigned kinematical groups. 
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Table 7. Contingency table between crossing groups and L, SR and M variability and M and C spectral types 





L-C 


SRb-C 


SRa-C 


M-C 


L-O 


SRb-0 


SRa-0 


M-0 


D(D) 


29 


25 


7 


9 


20 


29 


4 


31 


OD(D) 


1 


3 





1 


1 


4 


1 


9 


D(ODb) 


7 


3 


1 


4 


6 


33 


7 


28 


D(ODf) 


2 


4 





2 


42 


39 


3 


2 


OD(ODf) 


3 


5 





1 


34 


22 


2 





OD(ODb) 





3 








14 


28 


13 


54 


ED(ED) 


1 


1 





1 


4 


7 


3 


18 



Table 8. Attraction-repulsion indices between crossing groups and L, SR and M variability and M and C spectral 
types 







L-C 


SRb-C 


SRa-C 


M-C 


L-O 


SRb-0 


SRa-0 


M-0 


Diskl 


D(D) 


2.7 


2.3 


3.5 


2.0 


0.6 


0.7 


0.4 


0.9 


Diskl 


OD(D) 


0.7 


1.4 





1.6 





0.9 


0.9 


1.9 


Disk2 


D(ODb) 


0.3 


0.5 


0.8 


1.5 


0.3 


1.4 


1.4 


1.4 


Disk2 


D(ODf) 


0.3 


0.6 





0.7 


2.3 


1.5 


0.6 


0.1 


Old Disk 


OD(ODf) 


0.6 


1.1 





0.5 


2.5 


1.2 


0.4 





Old Disk 


OD(ODb) 





0.4 








0.6 


0.9 


2.2 


2.1 


Ext. Disk 


ED(ED) 


0.4 


0.4 





0.9 


0.6 


0.7 


1.6 


2.2 



(Mennessier et al., 1999) and with both observations 
and model predictions according to which dredge-up of 
carbon from core to the surface of AGB stars is more 
efficient in massive stars (Dopita et al.,1997). However, 
we must be ver y cau tions. Indeed, the biases described 
before (section |5. 2. 3 ) show the over-representativity of 
C-rich stars in the sample and the large number of 
missing LPVs. To which galactic population do the 
not-HIPPARCOS-observed LPVs belong? Are they O 
or C-rich? The ratio of carbon to oxygen-rich LPVs 
that goes up to 78% for the D(D) sample could be due 
to the sampling. Thus deduced implications for the de- 
pendence on the mass of the efficiency of the dredge-up 
are to be taken with caution. 

0-rich SRa's are distributed close to 0-rich Miras. 
0-rich SRb's are significantly present in the two disk 
population groups. The first one is composed of stars 
without a shell or with a thin one (cf sect.|4.3|). This 
probably corresponds to early AGB stars with a rela- 
tively high initial mass. The second one contains stars 
with a shell and corresponds to stars in the same stage 
as SRa and Miras. This agrees with results by Ker- 
schbaum and Hron (1992). 0-rich irregular LPVs are 
close to the SRb's lacking a shell and are probably 
early AGB stars. 

The difference between L irregular variable stars ac- 
cording to their spectral type is evident. C-rich L stars 
correspond to massive TP-AGB's with a thick enve- 
lope. On the contrary, 0-rich L stars are close to early 
AGB 0-rich SRb's but their initial mass range seems 
to be more extended to lower masses. 



One CH star (V Ari) is assigned to the Extended Disk 
population and it is the faintest star in all the luminosi- 
ties. This completely agrees with the usual hypothesis 
considering this type of star as an old giant star. Its 
amplitude of variability is less than 1 magnitude, a 
very small value even for a semi-regular. The C char- 
acter of this type of star is however difficult to explain. 
HIPPARCOS observations suggest that V Ari could be 
a suspected multiple star but this is not conclusive. 



6.5. Upper limit of the AGB 



In Section ^.2.2| we remarked that the brightest stars in 
the sample agree with the brightest luminosity for each 
group population (see figure |^) . Thus we can consider our 
sample as representative of the LPVs population as far as 
the brightest luminosities are concerned. Our calibrations 
show that the upper limit in K luminosity of the OD pop- 
ulation {{Kq — 3(Tif) = —8.1 mag.) is fainter than that 
of the D population {{Kq — Scix) — —9.4 mag.) as seen 
in Table ^j. This confirms the dependence of the upper 
limit of the AGB on Mms- Willson (1980) has described 
a schematic evolution on the AGB related to the mass-loss 
rate, its acceleration by the pulsations and probably the 
induced dust formation. She found a difference in solar lu- 
minosities of ^ 0.3L/ Lq where stars of solar abundance 
and Aims equal to 1.5 and 1 Mq leave the AGB. Our 
result is of the same order. 
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7. Conclusion 

Using available HIPPARCOS data wc apply the LM al- 
gorithm to improve the luminosity calibrations in visible, 
near-infrared and infrared wavelength ranges and to get 
information about the star and the circumstellar envelope. 

According to the galactic population - related to 
initial mass and mctallicity of the stars - and to the 
circumstellar envelope thickness and expansion, several 
groups of LPVs are obtained: bright (BD) and disk 
(diskl) galactic population with bright and expanding 
envelope, not so young and massive disk population 
(disk2) divided into 2 groups: one with thin envelope (f) 
and the other with a bright and expanding envelope (b). 
A similar separation according to envelope properties 
is found for the old disk (OD) population. At least 
some LPVs are found to belong to extended disk (ED) 
population. 

Our results deduced from kinematic properties con- 
firm that the AGB evolution depends on the initial mass 
of the progenitor in the main sequence. This agrees 
with the comparison of color-magnitude diagrams using 
our estimated individual luminosities with theoretical 
evolutionary tracks. According to the assigned galactic 
population wc can give ranges of age and of the lower 
limit main sequence mass for each star of our sample. 
The upper limit of the AGB also depends on Mms- The 
difference of the values found in K luminosity limits are 
consistent with Willson's schematic model related to the 
mass loss rate and its acceleration by the pulsations: 
"Stars evolve up the AGB with only moderate mass loss; 
at Te ~ 3000K Mira pulsation commences, driving the 
mass loss rate up by at least a factor 10". The induced 
dust formation is followed by the stabilization of the K 
luminosity after the carbon enrichment. 

The ultimate aim of this work is to estimate in- 
dividual K, 12 and 25 absolute magnitudes given, in 
the annex (available as an electronic table at CDS and 
in the ASTRID database). This allows us to study 
simultaneously the stellar properties and the behavior 
of the circumstellar envelope. The results recalled in the 
previous paragraph are obtained thanks to the estimated 
individual luminosities and they mainly concern proper- 
ties related to the assigned galactic populations. They 
will be systematically used in another paper (Mennessier 
et al., 2001) to study imphcations regarding the physics 
of LPVs, specifically the simultaneous stellar and cir- 
cumstellar evolutions along the Asymptotic Giant Branch. 
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